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1.  INTRODUCTION 


The  erosion  of  the  metallic  surface  on  the  inside  of  the  gun  tube  during  the  course 
of  firing  can  detrimentally  affect  the  performance  and  useful  life  of  the  gun  tube.  Gun 
tube  erosion  can  be  a  complex  phenomenon  because  there  are  many  possible  mechanisms 
that  contribute  to  the  erosion  of  the  inner  surface  of  the  gun  tube.  One  of  the  primary 
mechanisms  stems  from  the  heat  load  transferred  from  the  propellant  gases  to  the  gun  tube 
wall.  If  the  heat  load  is  large  enough,  melting  of  metal  at  the  inner  surface  of  the  gun 
tube  can  result.  The  high  velocity  propellent  gases  provide  a  means  of  removing  the  melted 
material  as  it  is  formed.  The  process  of  melting  and  subsequent  removal  of  the  melted 
material  is  referred  to  as  ablation.  The  ablation  process  as  a  potential  erosion  mechanism  is 
the  focus  of  the  current  numerical  study. 

Evidence  exists  from  previous  studies  that  ablation  is  a  potential  erosion  mechanism.  A 
recent  study  by  Bundy,  Gerber,  and  Bradley1  examined  the  single-shot  thermal  response  of 
the  M256  cannon  firing  the  M829A1  projectile.  Their  results  showed  that  the  chrome  plating 
on  the  inner  surface  of  the  gun  tube  provided  sufficient  protection  to  the  gun  tube  so  that 
melting  of  the  chrome  or  the  base  (underlying)  metal  would  not  occur  as  long  as  the  chrome 
plating  on  the  inner  surface  of  the  gun  tube  remained  intact.  Their  results  also  indicated 
that  if  the  chrome  layer  were  removed  (due  to  chipping  or  other  mechanical  means),  the  base 
metal  would  reach  its  melt  temperature  at  certain  locations  along  the  gun  tube.  Because 
their  numerical  approach  did  not  model  the  phase  change  of  the  metal,  no  predictions  of  the 
ablation  rate  could  be  obtained.  The  purpose  of  the  current  study  is  to  examine  the  ablation 
of  the  base  metal,  given  the  chrome  layer  has  been  removed  during  the  course  of  previous 
firings. 

In  order  to  examine  the  ablation  that  may  occur  during  heating  of  the  gun  tube,  a  one¬ 
dimensional  time-dependent  heat  conduction  model  has  been  developed  and  is  presented 
here.  The  code  uses  a  Crank-Nicolson  implicit  scheme  to  solve  the  governing  equation, 
which  is  cast  in  generalized  coordinate  form.  Prior  to  the  onset  of  melting,  heating  of  the 
inner  surface  of  the  gun  tube  is  modeled  using  a  convective  heating  boundary  condition. 
Once  the  melt  temperature  at  the  inner  surface  of  the  gun  tube  is  reached,  an  ablation 
boundary  condition  is  applied.  (Here,  ablation  occurs  when  the  solid  surface  reaches  the 
melt  temperature,  undergoes  phase  change  to  liquid,  and  the  liquid  phase  is  immediately 
removed  by  external  forces  such  as  shear  from  the  gun  gases.)  Since  the  model  is  a  one¬ 
dimensional  model,  only  radial  temperature  variations  are  considered.  Surface  heat  transfer 
coefficients  and  gas  temperatures  used  as  inputs  to  the  calculation  were  computed  externally 
from  the  heat  transfer  computation  using  an  interior  ballistics  code.  The  computed  heat 
transfer  coefficients  assume  that  the  inner  surface  of  the  gun  tube  is  smooth  and  do  not 
include  the  local  variations,  which  might  be  caused  by  erosion  cratering.  Further  details  on 
the  governing  equations  and  numerical  techniques  are  discussed  in  the  following  sections. 

The  numerical  technique  was  validated  with  previously  published  numerical  and  an¬ 
alytical  results  for  a  model  ablation  problem.  The  validated  technique  was  then  applied 
to  examine  the  thermal  performance  and  ablation  phenomenon  associated  with  the  M256 
cannon  firing  the  M829A1  kinetic  energy  projectile.  Predictions  of  the  in-depth  temperature 
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response  of  the  gun  tube  as  well  as  the  ablation  rate  were  made.  Details  of  the  validation 
and  application  are  given  in  the  results  section. 


2.  GOVERNING  EQUATIONS  AND  BOUNDARY 

CONDITIONS 


The  in-depth  temperature  response  of  the  unablated  (solid)  material  was  modeled  using 
the  one- dimensional  heat  conduction  equation  shown  in  equation  1. 


dT 

pCpdt  ~ 


r0  Qr 


(i) 


By  setting  f3  =  0  or  /3  =  1,  the  planar  or  axisymmetric  form  of  the  governing  equation 
can  be  obtained.  In  this  form  of  the  equation,  the  relevant  material  properties-the  density 
(p),  specific  heat  (cp),  and  the  conductivity  (k)- may  be  variable  (but  continuous).  Only 
constant  material  properties  were  considered  in  the  results  presented  here. 

At  the  inner  surface  of  the  gun  tube,  two  separate  boundary  conditions  were  applied 
depending  on  whether  or  not  melting  of  the  surface  material  was  occurring.  If  the  surface 
temperature  was  less  than  the  melt  temperature,  a  convective  heat  transfer  boundary  con¬ 
dition  (Newton’s  Law  of  Cooling)  was  applied.  This  boundary  condition  simply  states  that, 
at  the  inner  surface  of  the  gun  tube,  the  energy  loss  from  the  gun  gases  due  to  convection  is 
balanced  by  the  energy  gain  into  the  gun  tube  by  conduction. 


h(Tg  -  Twall)  =  -k^  (2) 

The  convective  heat  transfer  from  the  gun  gases  can  be  represented  by  a  convective 
heat  transfer  coefficient,  h,  and  a  driving  potential,  which  is  the  difference  between  the 
surface  temperature,  Twaii,  and  the  gas  temperature,  Tg.  Tg  (along  with  other  additional 
interior  ballistic  parameters)  was  obtained  using  the  the  NOVA  code2.  Using  results  from 
the  NOVA  code,  h  was  obtained  from  the  Veritay  code3' 4  using  the  correlation  of  Stratford 
and  Beavers5.  Computations  for  the  firing  of  the  M829A1  kinetic  energy  projectile  from  the 
M256  cannon  were  performed  prior  to  the  heat  conduction  computation.  Twaii  was  obtained 
by  coupling  the  boundary  condition  to  the  heat  conduction  equation. 

The  variations  of  the  convective  heat  transfer  coefficient  and  the  gas  temperature  with 
time  at  0.70  m  from  the  breech  are  shown  in  Figures  1  and  2,  respectively.  Both  the 
convective  heat  transfer  coefficient  and  the  gas  temperature  show  a  rapid  rise  after  shot 
passage  (~3  ms).  The  convective  heat  transfer  coefficient  reaches  a  maximum  about  2  ms 
after  shot  passage,  followed  by  a  rapid  decrease  as  the  firing  cycle  proceeds.  The  convective 
heat  transfer  coefficient  and  the  gas  temperature  vary  along  the  axial  length  of  the  gun  tube 
with  the  highest  levels  of  heat  transfer  occurring  near  the  breech. 

Once  the  inner  surface  of  the  gun  tube  reaches  the  melt  temperature,  melting  of  the 
surface  material  is  assumed  to  occur.  After  fully  melting  (complete  change  of  phase  from 
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solid  to  liquid),  the  liquid  metal  is  immediately  removed  by  the  shearing  action  of  the  gun 
gases.  Because  surface  material  is  being  removed  during  the  ablation  process,  the  surface 
location  (or  the  ablation  rate)  becomes  an  additional  unknown  in  the  problem.  During  the 
ablation  process,  two  boundary  conditions  are  applied  as  shown  in  equations  3  and  4. 


Twall 


Tmelt 


OT 

h(Tg  —  Twaii )  —  k-^ 


(3) 

(4) 


The  first  boundary  condition  states  that  the  temperature  at  the  surface  where  the  phase 
change  is  occurring  is  equal  to  the  melt  temperature.  The  melt  temperature  is  assumed 
to  be  a  known  material  property.  The  second  boundary  condition  is  obtained  from  an 
energy  balance  at  the  melt  surface.  In  addition  to  the  heat  transfer  due  to  conduction  and 
convection,  the  additional  energy  required  due  to  the  change  of  phase  from  solid  to  liquid 
(latent  heat  of  melting)  is  also  included.  This  boundary  condition  allows  the  ablation  rate  of 
ds 

the  solid  material,  — ,  to  be  computed,  given  that  the  latent  heat  of  melting,  L,  is  a  known 
at 


material  property. 


Several  milliseconds  after  shot  passage,  the  convective  heat  transfer  that  drives  the 
thermal  heating  of  the  gun  tube  decreases,  and  the  ablation  ceases.  The  convective  boundary 
condition  is  then  applied  as  before  the  melting  occurred  with  the  temperature  at  the  inner 
surface  computed  by  coupling  the  heat  conduction  equation  and  the  convective  boundary 
condition.  As  implemented  in  the  numerical  algorithm,  the  completion  of  the  melting  process 
is  detected  when  the  ablation  rate  changes  sign. 


To  compute  the  in-depth  temperature  response  of  the  gun  tube,  an  additional  boundary 
condition  is  required  at  the  outer  surface  of  the  gun  tube.  Here,  a  convective  boundary 
condition  can  be  applied: 

h osiT outer— wall  Too)  —  k  (5) 


where  Touter-waii  is  the  temperature  at  the  outer  surface  of  the  gun  tube,  hoo  is  the  convective 
heat  transfer  coefficient,  and  Too  is  the  ambient  temperature.  Representative  values  for  hoo 
and  can  be  found  in  Reference  1. 


Because  the  ablation  process  occurs  during  the  first  several  milliseconds  of  the  firing 
process,  the  heat  can  only  penetrate  a  fraction  of  the  distance  from  the  inner  surface  of  the 
gun  tube  to  the  outer  surface  of  the  gun  tube.  In  this  case,  it  may  not  be  necessary  to 
model  the  radial  temperature  response  of  the  entire  gun  tube.  By  locating  the  outer  edge  of 
the  computational  domain  far  enough  away  so  that  the  temperature  response  is  unchanged 
during  the  heating  process,  a  smaller  computational  domain  can  be  analyzed.  In  this  case, 
one  of  two  boundary  conditions  may  be  applied:  a  constant  temperature  boundary  condition 
or  an  adiabatic  wall  boundary  condition. 

The  appropriate  depth  of  the  computational  domain  can  be  estimated  using  a  depth- 
of-penetration  analysis.6’  7  The  depth  of  penetration,  6,  is  the  approximate  distance  that  the 
heat  would  penetrate  in  a  given  time,  t,  and  is  a  function  of  the  thermal  diffusivity  of  the 

material,  ct.  _ 

6  =  y/l2cd  (6) 
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Using  a  factor  of  safety  of  2  or  3  will  place  the  outer  boundary  of  the  computational  domain 
far  enough  from  the  inner  surface  of  the  gun  tube  so  that  the  temperature  at  the  outer 
boundary  remains  unaffected  by  the  heating  of  the  gun  tube  during  the  simulation. 


3.  COMPUTATIONAL  APPROACH 


The  governing  equations  and  boundary  conditions  are  solved  using  an  implicit  finite- 
difference  technique.  Prior  to  the  onset  of  melting,  the  governing  equations  and  boundary 
conditions  are  linear,  and  solutions  are  obtained  in  a  direct  (noniterative)  fashion.  During 
the  melting  process,  the  equations  become  nonlinear  because  the  dimensions  of  the  compu¬ 
tational  domain  are  coupled  with  the  ablation  rate.  An  iterative  approach  is  utilized  during 
melting  to  appropriately  address  the  nonlinearity. 

Because  the  boundary  of  the  computational  domain  moves  during  the  ablation  event,  a 
transformed  version  of  the  governing  equation  is  employed  in  the  computational  approach. 
This  allows  the  equations  to  be  solved  in  a  fixed  computational  space  even  though  the  physical 
boundary  is  moving.  A  generalized  transformation  between  the  computational  coordinate,  £, 
and  the  physical  coordinate,  r,  is  utilized.  The  transformed  equations  are  shown  as  follows: 


fdT  ,  dT , 
pCp{~dt+^}  ~ 

1  ,  d  f  Hut  dTs 

(7) 

dr 

6  = 

111  =  11 
dr 

(8) 

dt 

tr  = 

£  I  *-* 
III 

1¥|~ 

(9) 

dt 


In  this  form,  the  nonlinear  nature  of  the  governing  equation  produced  by  the  moving  bound¬ 
ary  is  evident  because  the  metric  terms,  £r  and  £<,  are  not  constant  and  are  dependent  on 
the  ablation  rate. 

The  equations  are  solved  using  a  Crank-Nicolson  finite-difference  technique.  The  finite- 
difference  form  of  the  governing  equations  is  shown  in  equation  10. 


JV+l  JV+l mN+l  _  N~NqrN  /mN+l  rpN+l\  / ' 

Pj-  p'  . ^  (  i+1  Ag"1  j  -  C1  -  *)6lf  (- 


jnJV  _  rpN  ' 


+0 


qy+1 

(rf+y 

t  !«■ 


f  j-iN+1 


tjpN+l y 


A£2 


/mN+l  mN+l\  ) 

-TO 


(rf+iykf+tJr  |?+4 


( 


J^N  _  rpN  > 


A£2  )} 

TfL  i 


A£2 


)} 


(10) 
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The  superscripts  N  and  N  +  1  represent  the  previous  and  current  time  steps,  while  the 
subscripts  j  —  1  and  j  +  1  denote  spacial  grid  node  locations  immediately  adjacent  to  node 
j.  Coefficients  that  are  to  be  evaluated  between  spacial  grid  nodes  j  —  1  and  j  or  j  and 
j  +  1  are  indicated  by  the  subscripts  j  —  \  and  j  +  |,  respectively.  The  parameter  0  allows 
for  switching  between  fully  explicit  and  fully  implicit  schemes.  In  the  computational  results 
presented  here,  the  second-order  accurate  Crank-Nicolson  scheme  (6  =  |)  was  utilized. 

Finite-difference  forms  of  the  boundary  conditions  were  also  used.  The  conduction 
terms  were  evaluated  using  a  first-order  one-sided  differencing  to  obtain  the  temperature 
gradient  at  the  inner  surface  of  the  gun  tube. 

When  melting  is  ongoing  and  the  grid  is  moving,  the  solution  is  iterated  at  each  time 
step  because  the  temperature,  ablation  rate,  grid  node  locations,  and  the  metric  terms  are 
coupled  in  a  nonlinear  fashion.  At  each  time  step,  the  location  of  the  melt  surface  at  the 
N  +  1  time  step  is  first  estimated  based  on  the  ablation  rate  at  the  previous  time  step.  The 
grid  node  locations  (which  are  dependent  on  the  melt  surface  location)  and  metric  terms  at 
the  N  +  1  time  step  axe  then  computed,  and  the  temperature  field  at  the  N  +  1  time  step 
can  then  be  updated  from  the  governing  equation  and  boundary  conditions.  Based  on  the 
predicted  temperature  field,  the  ablation  rate  at  the  N  +  1  time  step  is  determined  from  a 
finite- difference  form  of  equation  4.  Using  the  updated  ablation  rate,  a  new  iteration  loop 
can  be  started,  to  again  predict  the  temperature  field  and  ablation  rate  at  the  N  +  1  time 
step  until  both  the  temperature  field  and  the  ablation  rate  converge. 

To  resolve  the  large  temperature  gradients  at  the  inner  surface  of  the  gun  tube,  clus¬ 
tering  of  the  computational  grid  was  performed  using  an  exponential  stretching  function. 
The  resulting  grid  has  fine  spacing  at  the  inner  surface  of  the  gun  tube  and  progressively 
larger  spacing  away  from  the  wall.  The  grid  stretching  utilized  here  allowed  the  grid  spacing 
between  adjacent  nodes  to  increase  by  no  more  than  10%  because  highly  stretched  grids  are 
known  to  reduce  numerical  accuracy. 

The  depth  of  penetration  analysis,  discussed  earlier,  can  also  be  applied  to  select  the 
spacing  of  the  computational  grid  at  the  inner  wall  of  the  gun  tube.  After  selecting  the 
time-step  for  the  integration  procedure,  the  grid  spacing  at  the  inner  wall  of  the  gun  tube 
can  be  estimated  by  computing  the  depth  of  penetration  after  one  time  step.  Numerical 
experiments  indicated  that  the  grid  spacing  at  the  wall  should  be  less  that  10%  of  the  depth 
of  penetration  after  one  time  step  to  avoid  inaccuracies  in  the  predicted  temperature  profile 
and  ablation  rate. 


4.  RESULTS 


The  numerical  method  has  been  validated  using  previous  numerical  and  analytical  ap¬ 
proaches  for  model  problems  involving  phase  change.  These  results  are  presented  in  the  next 
section.  After  validating  the  numerical  approaches,  application  of  the  technique  was  made 
to  examine  the  ablation  phenomenon  in  the  M256  cannon  firing  M829A1  kinetic  energy 
projectiles.  Presentation  of  these  results  follows  the  validation  results. 
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4.1  Validation:  Ablation  of  a  Semi-Infinite  Solid 


Over  the  past  several  decades,  numerous  studies  examining  the  phase  change  process 
have  been  made.  Results  from  two  of  these  studies  have  been  used  to  validate  the  currently 
developed  numerical  method.  Landau8  has  made  one-dimensional  time-dependent  numerical 
predictions  of  a  melting  solid  using  a  technique  similar  to  that  utilized  here.  Landau  consid¬ 
ered  the  case  of  a  semi-infinite  solid  subjected  to  a  constant  heat  flux.  Goodman9  considered 
similar  problems  using  a  heat-balance  integral  approach.  The  approach  of  Goodman  utilizes 
an  assumed  form  of  the  temperature  profile  to  analytically  determine  the  heat  conduction 
and  ablation  process.  While  exact  solutions  axe  not  typically  obtained,  the  results  axe  rea¬ 
sonably  accurate  and  axe  of  a  simple  form.  One  of  the  cases  addressed  by  Goodman  was  the 
melting  of  a  semi-infinite  solid  under  constant  heat  flux.  In  both  cases,  an  ablative  boundary 
condition  was  utilized;  that  is,  the  liquid  phase  was  immediately  removed  following  melting. 

The  problem  of  the  heating  on  a  semi-infinite  solid  subjected  to  constant  heat  flux  was 
also  addressed  using  the  current  numerical  method.  A  finite  computational  domain  was 
actually  utilized  in  the  computation,  though  the  unheated  boundary  was  placed  far  enough 
from  the  heated  boundary  so  that  the  temperature  near  the  unheated  boundary  did  not 
change  during  the  computation.  The  solid  initially  had  a  uniform  temperature  distribution 
of  T0  and  a  specified  melt  temperature,  Tmeit •  The  results  are  scaled  by  the  time  to  melt, 
tmeit,  which  is  a  function  of  applied  heat  flux.  The  instantaneous  location  of  the  solid  surface, 
s(t),  relative  to  the  fixed  coordinate,  r,  is  shown  schematically  in  Figure  3. 

Figure  4  shows  the  predicted  in-depth  temperature  profile  at  four  different  times  during 
the  ablation  process.  The  results  were  obtained  using  material  properties  corresponding  to 
gun  steel  (Table  1).  In  addition  to  the  numerical  results,  the  analytical  results  obtained 
using  the  heat-balance  integral  approach  of  Goodman  are  also  shown  at  the  onset  of  melting 
and  during  steady-state  melting.  The  numerical  solution  and  the  analytical  results  differ  by 
less  than  0.2%  at  the  onset  of  melting.  (Though  not  shown  both  results  also  show  similar 
agreement  with  the  exact  analytical  solution.)  As  the  melting  progresses,  the  tempera¬ 
ture  gradient  at  the  surface  is  reduced  due  to  the  additional  energy  required  to  melt  the 
solid.  After  approximately  eight  times  the  time  required  to  reach  the  melt  temperature, 
the  numerical  results  show  that  the  temperature  profile  has  nearly  reached  the  steady-state 
temperature  profile.  The  differences  between  the  analytical  and  numerical  steady-state  tem¬ 
perature  profile  are  less  than  2%  of  the  difference  between  the  melt  temperature  and  the 
initial  temperature.  Exact  agreement  between  the  two  results  is  not  expected  because  the 
assumed  form  of  the  temperature  profile  in  the  heat-balance  integral  approach  is  not  an 
exact  solution. 

A  comparison  of  the  predicted  ablation  rate  obtained  using  the  current  technique  with 
the  results  of  Landau  is  shown  in  Figure  5.  The  computed  ablation  rates  are  scaled  by 
the  steady-state  ablation  rate,  which  can  be  calculated  analytically.  The  current  results 
are  generally  -within  l%-2%  of  the  results  of  Landau  for  the  three  different  sets  of  material 
properties,  including  properties  close  to  gun  steel.  The  largest  difference  between  the  two 
sets  of  results  is  4%.  (Differences  between  the  results  of  Landau  and  the  current  result  may 
be  more  related  to  problems  in  extracting  the  data  from  the  original  published  graphs  of 
Landau  than  to  numerical  accuracy.) 
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The  results  show  that  the  ablation  rate  approaches  the  steady-state  ablation  rate  in  an 
asymptotic  manner  with  materials  having  a  larger  value  of  m  approaching  the  steady-state 
ablation  rate  at  earlier  nondimensional  times.  The  predictions  of  the  steady-state  ablation 
rate  were  obtained  by  running  the  computation  for  extended  time  (t  >  15 tmeit  —  200 tmeit 
depending  on  the  material  property  m).  The  predicted  values  were  within  0.1%  of  the  exact 
analytical  results  presented  by  Landau. 


4.2  Application  to  the  M256  Cannon  Firing  an  M829A1  Projectile 

Application  of  the  numerical  technique  for  predicting  gun  tube  erosion  was  made  to  the 
M256  gun  tube  firing  an  M829A1  projectile.  The  gun  tube  is  fabricated  with  AISI 4340  steel 
with  the  properties  shown  in  Table  1.  The  inner  surface  of  the  gun  tube  is  coated  with  a 
thin  chrome  layer  which  protects  the  gun  tube.  The  chrome  has  a  higher  melt  temperature 
and  hardness  than  the  gun  steel. 

Preliminary  computations  demonstrated  that  with  the  chrome  layer  intact,  both  the 
chrome  layer  and  the  gun  steel  would  not  melt  (or  ablate)  due  to  convective  heating  during 
a  single  shot  firing  cycle.  However,  if  the  chrome  layer  were  removed,  the  inner  surface  of 
the  gun  tube  could  reach  its  melt  temperature  during  the  firing  cycle  at  locations  near  the 
breech  of  the  gun.  This  behavior  has  been  noted  previously1.  It  has  been  noted  that  the 
chipping  of  the  chrome  layer  can  occur  after  several  hundred  firings.  The  results  presented 
here  address  the  erosion  of  the  gun  tube  following  the  removal  of  the  chrome  layer. 

Computations  of  the  thermal  response  of  the  gun  tube  were  made  for  a  single-shot 
firing  scenario  with  the  initial  gun  tube  temperature  assumed  to  be  294  K.  The  propelling 
charge  temperature  was  294  K  corresponding  to  ambient  conditions.  The  thermal  response 
of  the  gun  tube  was  computed  at  several  axial  locations  along  the  gun  tube  to  demonstrate 
variation  in  the  thermal  response  along  the  length  of  the  gun  tube. 

The  temperature  response  of  the  gun  tube  on  its  inner  surface  as  a  function  of  time  is 
shown  in  Figure  6.  Results  at  four  different  axial  stations  along  the  gun  tube  are  shown. 
At  0.70  m  and  0.76  m  from  the  breech,  the  inner  surface  of  the  gun  tube  reaches  the  melt 
temperature,  and  some  ablation  of  the  inner  surface  occurs.  Further  down  the  bore  of  the 
gun  tube,  the  convective  heat  transfer  drops  to  a  level  where  melting  of  the  gun  steel  does 
not  occur. 

Since  the  ablation  occurs  over  such  a  short  duration,  it  appears  to  have  little  effect  on 
the  temperature  response  of  the  gun  tube  except  during  the  time  period  when  the  ablation 
occurs.  Figure  7  shows  the  temperature  response  of  the  gun  tube  at  0.70  m  from  the  breech. 
Two  computational  results  are  shown:  one  with  the  ablation  model  and  one  with  no  ablation 
model.  The  results  differ  by  less  than  20  K  immediately  after  the  ablation  is  complete,  with 
the  difference  decreasing  to  less  than  5  K  later  in  the  firing  cycle. 

In  other  applications,  such  as  thermal  protection  for  spacecraft  re-entry,  the  ablation 
process  can  provide  a  measure  of  thermal  protection,  in  that  a  portion  of  the  heat  input 
into  the  solid  is  expended  in  the  phase  change  of  the  material.  However,  due  to  the  short 
duration  of  the  ablation  process  in  this  application,  the  ablation  does  not  have  much  effect 
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on  the  total  heat  input  into  the  gun  during  a  single-shot  firing.  For  example,  at  0.70  m 
from  the  breech,  the  difference  in  the  total  heat  input  to  the  gun  tube  during  firing  with  or 
without  the  ablation  model  is  less  than  1%. 

The  ablation  rate  of  the  gun  tube  as  a  function  of  time  at  two  locations  near  the  breech  is 
shown  in  Figure  8.  The  location  0.70  m  from  the  breech  shows  a  much  higher  rate  of  ablation 
than  at  0.76  m  from  the  breech.  As  shown  in  the  previous  figure,  this  location  reaches  the 
melt  temperature  first  and  remains  at  the  melt  temperature  for  a  longer  duration.  The  total 
ablation  for  a  single  firing  was  1.97  fim  and  0.31  (xm  at  0.70  m  and  0.76  m  from  the  breech, 
respectively.  Over  400  firings,  this  ablation  rate  would  translate  to  a  total  eroded  depth  of 
0.8  mm  and  0.12  mm  at  these  two  stations.  The  computed  erosion  depths  are  apparently 
the  right  order  of  magnitude  for  normal  observed  erosion1 . 

The  in-depth  temperature  profile  during  the  ablation  event,  shown  in  Figure  9,  demon¬ 
strates  the  limited  depth  of  heat  penetration  since  the  beginning  of  the  firing  cycle.  The 
results  justify  the  use  of  a  one- dimensional  approach  for  examining  the  ablation  event  because 
the  temperature  profile  should  be  influenced  by  the  conditions  immediately  surrounding  a 
given  location.  The  results  also  imply  that  some  types  of  active  measures  for  removing  the 
heat  from  the  gun  tube  may  be  ineffective  for  preventing  erosion  unless  they  can  affect  the 
conditions  very  close  to  the  gun  tube  wall. 

Clearly,  the  use  of  multidimensional  approaches  would  be  required  if  the  geometric 
effects  resulting  from  gun  tube  rifling  or  erosion  channels  axe  to  be  investigated.  The  devel¬ 
opment  of  such  approaches  is  seen  as  a  desirable  research  objective. 

Additional  computational  results  were  generated  to  quantify  the  sensitivity  of  the  ab¬ 
lation  process  to  the  convective  heating.  The  results  showed  that  a  4%  reduction  in  the 
convective  heat  transfer  coefficient  or  a  2%  reduction  in  the  gun  gas  temperature  is  enough 
to  eliminate  the  ablation  at  0.70  m  from  the  breech.  This  is  one  indicator  of  how  sensitive 
the  ablation  event  might  be  to  effects  that  might  augment  or  retard  the  heating  at  the  in¬ 
bore  surface.  It  also  serves  as  a  caution  to  the  analyst  who  is  modeling  erosion  because  even 
though  the  bulk  of  the  input  heat  may  be  provided  by  the  convective  heating,  additional 
minor  sources  of  heating  may  yield  enough  additional  heat  to  contribute  significantly  to  the 
total  amount  of  erosion. 


5.  CONCLUSION 


A  numerical  method  for  estimating  the  erosion  that  occurs  within  gun  tubes  due  to  high 
convective  in-bore  heating  has  been  developed.  The  method  has  been  successfully  applied  to 
estimate  the  erosion  that  may  occur  during  the  firing  of  M829A1  kinetic  energy  projectiles 
fired  from  the  M256  cannon,  the  main  gun  on  the  Ml  tank. 

Numerical  predictions  have  demonstrated  that,  with  the  chrome  layer  intact,  melting 
or  ablation  of  the  chrome  layer  or  the  base  gun  steel  is  not  expected  to  occur.  However,  the 
predictions  show  that  if  the  chrome  layer  is  removed  (which  has  been  observed  after  several 
hundred  firings),  erosion  can  occur  within  a  region  between  the  breech  and  0.75  m  down- 
bore  of  the  breech.  The  predictions  show  that  the  erosion  occurs  over  a  very  short  duration 
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shortly  after  shot  passage.  The  predicted  levels  of  erosion  are  consistent  with  experimental 
observations. 

During  the  course  of  the  study,  one  significant  technical  issue  emerged:  the  accurate 
prediction  of  erosion  rates  may  be  a  challenging  endeavor  because  the  erosion  occurs  over 
a  very  short  duration  and  the  rate  is  quite  sensitive  to  a  number  of  factors.  In  this  study, 
care  was  taken  to  ensure  grid-independent  solutions.  Inaccuracies  were  demonstrated  when 
inappropriate  selection  of  grid  parameters  such  as  grid  stretching  and  spacing  was  made. 
Other  issues  such  as  the  accuracy  of  the  predicted  surface  heat  transfer  were  not  investigated 
in  this  study  but  potentially  have  a  significant  effect  on  the  predicted  erosion  rate.  (Many 
other  studies  have  addressed  this  issue,  but  typically  these  studies  have  focused  on  the  total 
heat  transfer  to  the  gun  tube,  rather  than  the  peak  heat  transfer  that  really  drives  the 
erosion.)  While  these  may  be  issues  for  the  physical  modeler,  this  is  actually  beneficial  for 
the  ballistician  because  it  means  that  subtle  changes  in  the  design  can  dramatically  change 
the  erosion  characteristics  of  the  gun  tube/projectile  system.  The  development  of  analysis 
tools  for  gun  tube  erosion  can  provide  the  designer  with  important  information  as  to  whether 
design  changes  improve  or  degrade  the  erosion  characteristics  of  a  particular  design. 
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Table  1.  Thermal  Properties  of  Gun  Steel 


Gun  Steel 

Chrome 

Cp 

469.05  J/kg-K 

505.3  J/kg-K 

k 

38.07  J/(m-s-K) 

83.75  J/(m-s-K) 

p 

7827.  kg/m3 

7191.  kg/m3 

L 

270.  kJ/kg 

- 

Tmelt 

1780.  K 

2130.  K 

10 


Figure  1.  Predicted  convective  heat  transfer  coefficient  versus  time,  0.70  m  from  the  breech, 
M256  cannon  firing  an  M829A1  projectile. 


Figure  2.  Predicted  in-bore  gas  temperature  versus  time,  0.70  m  from  the  breech,  M256 
cannon  firing  the  M829A1  projectile. 
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Solid  Surface  Location 


Figure  3.  Schematic  for  ablation  of  a  semi-infinite  solid. 
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Figure  4.  In-depth  temperature  profile  during  ablation  process  for  a  semi-infinite  slab 
subject  to  constant  rate  heat  flux. 
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Figure  5.  Fractional  ablation  rate  versus  time  for  a  semi-infinite  slab  subject  to  constant 
rate  heat  flux. 


Figure  6.  In-bore  surface  temperature  of  gun  tube  during  firing  cycle  at  various  axial 
locations. 
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Figure  7.  In-bore  surface  temperature  of  gun  tube  during  firing  cycle  at  0.70  m  from  the 
breech,  with  and  without  ablation  model. 
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Figure  8.  Ablation  rate  of  in-bore  surface  of  gun  tube  during  firing  cycle  at  two  axial 
locations. 
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Figure  9.  In-depth  temperature  profile  of  gun  tube  0.70  m  from  breech  during  ablation 
process  (5  ms). 
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LIST  OF  SYMBOLS 


Cp 

ds 

dt 

h 

h<x> 

k 

L 

m 


s 

t 

tmelt 

T 

Too 

Ta 

Tmelt 

T0 

TyjO.ll 

Touter—wall 


specific  heat 
ablation  rate 

convective  heat  transfer  coefficient 

convective  heat  transfer  coefficient  on  exterior  of  gun  tube 

conductivity- 

latent  heat  of  melting 

heat  capacity  to  latent  heat  ratio 

radial  coordinate 

instantaneous  location  of  solid  surface 
time 

time  to  melt 

temperature 

ambient  temperature 

gas  temperature 

melt  temperature 

initial  temperature 

interior  wall  or  surface  temperature 

external  wall  or  surface  temperature 


Greek  Symbols 


a 

fi 

6 

6 

t 

P 


thermal  diffusivity 

switch  for  planar  (f3  =  0)  or  axisymmetric  (/?  =  1)  form  of  equations 

depth  of  heat  penetration 

temporal  differencing  parameter 

computational  coordinate 

density 
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CDR  US  ARMY  ARDEC 
ATTN  AMSTA  AR  AEE  B 
D  DOWNS 
S  EINSTEIN 
S  WESTLEY 
S  BERNSTEIN 
J  RUTKOWSKI 
B  BRODMAN 
P  O’REILLY 
R  CIRINCIONE 
PHUI 
J  O’REILLY 

PCNTY  ARSNLNJ  07806-5000 

CDR  US  ARMY  ARDEC 

ATTN  AMSTA  AR  AEE  WW 

MMEZGER 

J  PINTO 

D  WIEGAND 

PLU 

CHU 

PCNTY  ARSNLNJ  07806-5000 

CDR  US  ARMY  ARDEC 
ATTN  AMSTA  AR  AES 
s  KAPLOwrrz 

PCNTY  ARSNLNJ  07806-5000 

CDR  US  ARMY  ARDEC 
ATTN  AMSTA  AR  HFM 
EBARRIERES 

PCNTY  ARSNLNJ  07806-5000 

CDR  US  ARMY  ARDEC 
ATTN  AMSTA  AR  FSA  T 
M  SALSBURY 

PCNTY  ARSNLNJ  07806-5000 

CDR  US  ARMY  ARDEC 

ATTN  AMSTA  AR  FSA  F  LTC  R  RIDDLE 

PCNTY  ARSNLNJ  07806-5000 

CDR  US  ARMY  ARDEC 

ATTN  AMSTA  AR  FSC  G  FERDINAND 

PCNTY  ARSNLNJ  07806-5000 

COMMANDER 

ATTN  AMSTA  AR  FS  T  GORA 
PCNTY  ARSNLNJ  07806-5000 


1  CDR  US  ARMY  ARDEC 

ATTN  AMSTA  AR  FS  DH 
JFENECK 

PCNTY  ARSNLNJ  07806-5000 

3  CDR  US  ARMY  ARDEC 

ATTN  AMSTA  AR  FSS  A 
R  KOPMANN 
B  MACHEK 
LPINDER 

PCNTY  ARSNLNJ  07806-5000 

1  CDR  US  ARMY  ARDEC 
ATTN  AMSTA  AR  FSN  N 
K CHUNG 

PCNTY  ARSNLNJ  07806-5000 

2  DIR  BENET  WEAPONS  LABS 
ATTN  AMSTA  AR  CCB  RA 
GP  O’HARA 

G  APFLEGL 

WATERVLIET  NY  12189-4050 

1  DIR  BENET  WEAPONS  LABS 

ATTN  AMSTA  AR  CCB  RT 
S  SOPOK 

WATERVLIET  NY  12189-4050 

1  DIR  BENET  WEAPONS  LABS 
ATTN  AMSTA  AR  CCB  S 
FHEISER 

WATERVLIET  NY  12189-4050 

2  CDR  US  ARMY  RSRCH  OFC 
ATTN  TECHNICAL  LIBRARY 
D  MANN 

POBOX  12211 

RSCHTRIPKNC  27709-2211 

1  CDR  USACECOM 

R&D  TECHNICAL  LIBRARY 
ATTN  ASQNC  ELC  IS  L  R 
MYER  CENTER 

FT  MONMOUTH  NJ  07703-5301 

1  CMDT  USA  AVIATION  SCHL 

ATTN  AVIATION  AGENCY 
FT  RUCKER  AL  36360 

1  PM  US  TANK  AUTOMOTIVE  CMD 

ATTN  AMCPM  ABMS  T  DEAN 
WARREN  MI  48092-2498 
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1  PM  US  TANK  AUTOMOTIVE  CMD 

FIGHTING  VEHICLE  SYSTEMS 
ATTN  SFAE  ASM  BV 
WARREN  MI  48397-5000 

1  PM  ABRAMS  TANK  SYSTEM 

ATTN  SFAE  ASM  AB 
WARREN  MI  48397-5000 

1  DIR  HQ  TRAC  RPD 

ATTN  ATCD  MA 
FT  MONROE  VA  23651-5143 

1  COMMANDER 

US  ARMY  BELVOIR  R&D  CTR 
ATTN  STRBE  WC 
FT  BELVOIR  VA  22060-5006 

1  DIRECTOR 

US  ARMY  TRAC  FT  LEE 
ATTN  ATRC  L  MR  CAMERON 
FT  LEE  VA  23801-6140 

1  COMMANDANT 

US  ARMY  CMD  &  GEN 
STAFF  COLLEGE 
FT  LEAVENWORTH  KS  66027 

1  COMMANDANT 

USA  SPECIAL  WARFARE  SCHL 
ATTN  REV  AND  TRNG  LIT  DIV 
FT  BRAGG  NC  28307 

1  COMMANDER 

RADFORD  ARMY  AMMUNITION  PLANT 
ATTN  SMCAR  QA  HI  LIB 
RADFORD  VA  24141-0298 

1  COMMANDER 

US  ARMY  NGIC 
ATTN  AMXST  MC  3 
220  SEVENTH  STREET  NE 
CHRLTTESVLLE  VA  22901-5396 

1  COMMANDANT 

USAFAC&S 

ATTN  ATSF  CD  COL  T  STRICKLIN 
FT  SILL  OK  73503-5600 


1  COMMANDANT 
USAFC&S 
ATSF  CN  P  GROSS 
FT  SILL  OK  73503-5600 

1  CMDT  USA  ARMOR  SCHL 
ARMOR  AGENCY 

ATTN  ATZK  CD  MS 
M  FALKOVITCH 
FT  KNOX  KY  40121-5215 

2  CDR  NAVAL  SEA  SYSTEMS  CMD 

ATTN  SEA  62R 

SEA  64 

WASH  DC  20362-5101 

1  CDR  NAVAL  AIR  SYSTEMS  CMD 

ATTN  AIR  954  TECH  LIBRARY 
WASH  DC  20360 

4  CDR  NAVAL  RSRCH  LAB 

ATTN  TECHNICAL  LIBRARY 
CODE  4410 
K  KAILASANATE 
J  BORIS 
E  ORAN 

WASH  DC  20375-5000 

1  OFFICE  OF  NAVAL  RSRCH 

ATTN  CODE  473  R  S  MILLER 
800  N  QUINCY  STREET 
ARLINGTON  VA  22217-9999 

1  OFFICE  OF  NAVAL  TECHLGY 

ATTN  ONT  213  D  SIEGEL 
800  N  QUINCY  ST 
ARLINGTON  VA  22217-5000 

1  CDR  NAVAL  SURFACE  WARFARE 

CTR 

ATTN  CODE  730 

SILVER  SPRING  MD  20903-5000 

1  CDR  NAVAL  SURFACE  WARFARE 

CTR 

ATTN  CODE  R  13  R  BERNECKER 
SILVER  SPRING  MD  20903-5000 
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7  CDR  NAVAL  SURFACE  WARFARE 

CIR 

ATTN  T  C  SMITH 
K  RICE 
S  MITCHELL 
S PETERS 
J  CONSAGA 
C  GOTZMER 
TECHNICAL  LIBRARY 
INDIAN  HEAD  MD  20640-5000 

1  CDR  NAVAL  SURFACE  WARFARE 

cm 

ATIN  CODE  G30  GUNS  & 
MUNITIONS  DIV 
DAHLGRENVA  22448-5000 

1  CDR  NAVAL  SURFACE  WARFARE 

cm 

ATTN  CODE  G32 
GUNS  SYSTEMS  DIV 
DAHLGRENVA  22448-5000 

1  CDR  NAVAL  SURFACE  WARFARE 

cm 

ATTN  CODE  G33  T  DORAN 
DAHLGRENVA  22448-5000 

1  CDR  NAVAL  SURFACE  WARFARE 

cm 

ATTN  CODE  E23 
TECHNICAL  LIBRARY 
DAHLGRENVA  22448-5000 

2  CDR  NAVAL  AIR  WARFARE  CTR 
ATTN  CODE  388 

CF  PRICE 
T BOGGS 

CHINA  LAKE  CA  93555-6001 

2  CDR  NAVAL  AIR  WARFARE  CTR 

ATTN  CODE  3895 
T  PARR 
R  DERR 

CHINA  LAKE  CA  93555-6001 

1  CDR  NAVAL  AIR  WARFARE  CTR 
INFORMATION  SCIENCE  DIV 
CHINA  LAKE  CA  93555-6001 


1  COMMANDING  OFFICER 

ATIN  CODE  5B331  TECH  LIB 
NAVAL  UNDERWATER  SYS  CTR 
NEWPORT  RI  02840 

1  AFOSR  NA 

ATIN  J  TISHKOFF 
BOLLING  AFB  DC  20332-6448 

1  OLAC  PL  TSTL 

ATTN  D  SHIPLETT 
EDWARDS  AFB  CA  93523-5000 

3  ALLSCF 

ATIN  J  LEVINE 
L  QUINN 
T  EDWARDS 

EDWARDS  AFB  CA  93523-5000 

1  WLMNAA 

ATIN  B  SIMPSON 
EGLINAFBFL  32542-5434 

1  WLMNME 

ENERGETIC  MATERIALS  BR 
2306  PERIMETER  RD 
STE  9 

EGLINAFBFL  32542-5910 

1  WLMNSH 

ATTN  R  DRABCZUK 
EGLINAFBFL  32542-5434 

2  NASA  LANGLEY  RSRCH  CIR 
ATIN  M  S  408 

W  SCALLION 
D  WTTCOFSKI 
HAMPTON  VA  23605 

1  CENTRAL  INTELLIGENCE  AGENCY 

OFC  OF  THE  CENTRAL  REFERENCES 
DISSEMINATION  BRANCH 
ROOM  GE  47  HQS 
WASHINGTON  DC  20502 

1  CENTRAL  INTELLIGENCE  AGENCY 

ATIN  J  BACKOFEN 
NHB  ROOM  5N01 
WASHINGTON  DC  20505 
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1  SDIO  TNI 

ATTN  L  H  CAVENY 
PENTAGON 

WASHINGTON  DC  20301-7100 

1  SDIO  DA 
ATTN  E  GERRY 
PENTAGON 

WASHINGTON  DC  21301-7100 

2  HQ  DNA 
ATTN  D  LEWIS 
A  FAHEY 

6801  TELEGRAPH  RD 
ALEXANDRIA  VA  22310-3398 

1  DIR  SANDIA  NATL  LABS 

ATTN  M  BAER 
DEPARTMENT  1512 
PO  BOX  5800 

ALBUQUERQUE  NM  87185 

1  DIR  SANDIA  NATL  LABS 

ATTN  R  CARLING 
COMBUSTION  RSRCH  FACILITY 
LIVERMORE  CA  94551-0469 

1  DIR  SANDIA  NATL  LABS 
ATTN  8741  G  A  BENEDOTT 
PO  BOX  969 

LIVERMORE  CA  94551-0969 

2  DIR  LLNL 
ATTN  L  355 

A  BUCKINGHAM 
M  FINGER 
PO  BOX  808 

LIVERMORE  CA  94550-0622 

1  DIRLANL 

ATTN  T3  D  BUTLER 

PO  BOX  1663 

LOS  ALAMOS  NM  87544 

1  DIRLANL 

ATTN  M  DIVISION  B  CRAIG 

PO  BOX  1663 

LOS  ALAMOS  NM  87544 


2  BATTELLE 

ATTN  TWSTTAC 
V  LEVIN 

505  KING  AVENUE 
COLUMBUS  OH  43201-2693 

1  BATTELLE  PNL 

ATTN  MR  MARK  GARNICH 
PO  BOX  999 
RICHLAND  WA  99352 

1  THE  UNIV  OF  AUSTIN  TEXAS 
INST  FOR  ADVANCED  TECHLGY 
ATTN  T  M  KIEHNE 

4030  2  W  BRAKER  LANE 
AUSTIN  IX  78759-5329 

2  CPIAJHU 

ATTN  H  J  HOFFMAN 
T  CHRISTIAN 

10630  LITTLE  PATUXENT  PKWY 
SUITE  202 

COLUMBIA  MD  21044-3200 

1  AFELM  THE  RAND  CORP 

ATTN  LIBRARY  D 
1700  MAIN  STREET 
SANTA  MONICA  CA  90401-3297 

1  BRIGHAM  YOUNG  UNIV 

ATTN  M  BECKSTEAD 
DEPT  OF  CHEMICAL  ENGRG 
PROVO  UT  84601 

1  CALIF  INSTITUTE  OF  TECHLGY 

ATTN  L  D  STRAND  MS  125  224 
JET  PROPULSION  LAB 
4800  OAK  GROVE  DRIVE 
PASADENA  CA  91109 

1  CALIF  INSTITUTE  OF  TECHLGY 

ATTN  F  E  C  CULICK 
204  KARMAN  LAB 
MAIN  STOP  301  46 
1201  E  CALIFORNIA  SIR 
PASADENA  CA  91109 
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3  GEORGIA  INSTITUTE  OF  TECHLGY 
SCHOOL  OF  AEROSPACE  ENGRG 
ATTN  B  T  ZIM 

E  PRICE 
WC  STRAHLE 
ATLANTA  GA  30332 

2  UNTV  OF  ILLINOIS 

DEPT  OF  MECH  INDUSTRY  ENGRG 
ATTN  H  KRIER 
R  BEDDINI 

144  MEB  1206  N  GREEN  ST 
URBANAIL  61801-2978 

1  UNTV  OF  MASSACHUSETTS 

DEPT  OF  MECHANICAL  ENGRG 
ATTN  K  JAKUS 
AMHERST  MA  01002-0014 

1  UNTV  OF  MINNESOTA 

DEPT  OF  MECHANICAL  ENGRG 
ATTN  E  FLETCHER 
MINNEAPOLIS  MN  55414-3368 

4  PENNSYLVANIA  STATE  UNIV 
DEPT  OF  MECHANICAL  ENGRG 
ATTN  V  YANG 

K  KUO 
C  MERKLE 
G  SETTLES 

UNIVERSITY  PARK  PA  16802-7501 

1  RENSSELAER  POLYTECH  INST 
DEPT  OF  MATHEMATICS 
TROY  NY  12181 

1  STEVENS  INSTIT  UTE  OF  TECHLGY 

ATTN  R  MCALEVY  HI 
DAVIDSON  LABORATORY 
CASTLE  POINT  STATION 
HOBOKEN  NJ  07030-5907 

1  RUTGERS  UNIVERSITY 

DEPT  OF  MECH  AND  AERO  ENGRG 
ATTN  S  TEMKIN 
UNIVERSITY  HEIGHTS  CAMPUS 
NEW  BRUNSWICK  NJ  08903 

1  UNIV  OF  UTAH 

DEPT  OF  CHEMICAL  ENGRG 
ATTN  A  BAER 

SALT  LAKE  CITY  UT  84112-1194 


1  WASHINGTON  STATE  UNIV 

DEPT  OF  MECHANICAL  ENGRG 
ATTN  C  T  CROWE 
PULLMAN  WA  99163-5201 

1  ARROW  TECHLGY  ASSOC  INC 
ATTN  W  HATHAWAY 

PO  BOX  4218 

SOUTH  BURLINGTON  VT 

05401-0042 

2  AAI  CORPORATION 
ATTN  JFRANKLE 
D  CLEVELAND 

PO  BOX  126 

HUNT  VALLEY  MD  21030-0126 

8  ALLIANT  TECHSYSTEMS  INC 

ATTN  R  E  TOMPKINS 
J KENNEDY 
J  BODE 
C  CANDLAND 
L  OSGOOD 
RBURETTA 
R BECKER 
M  SWENSON 
600  SECOND  ST  NE 
HOPKINS  MN  55343 

1  ELI  FREEDMAN  AND  ASSOCIATES 

ATTN  E  FREEDMAN 
2411  DIANA  RD 
BALTIMORE  MD  21209-1525 

1  GENERAL  APPLIED  SCIENCES  LAB 

ATTN  J  ERDOS 
77  RAYNOR  AVE 
RONKONKAMANY  11779-6649 

1  GENERAL  ELECTRIC  COMPANY 

TACTICAL  SYSTEM  DEPT 
ATTN  JMANDZY 
100  PLASTICS  AVE 
PITTSFIELD  MA  01201-3698 

1  HTRI 

ATTN  M  J  KLEIN 
10  W  35TH  STREET 
CHICAGO  EL  60616-3799 
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HERCULES  INC 
ATTN  L  GIZZI 
D  A  WORRELL 
W  J  WORRELL 
C CHANDLER 

RADFORD  ARMY  AMMO  PLANT 
RADFORD  VA  24141-0299 

HERCULES  INC 
ATTN  WILLIAM  B  WALKUP 
THOMAS  F  FARABAUGH 
ALLEGHENY  BALLISTICS  LAB 
PO  BOX  210 

ROCKET  CENTER  WV  26726 


1  OLIN  ORDNANCE 

ATTN  H  A  MCELROY 
10101  9TH  STREET  NORTH 
ST  PETERSBURG  FL  33716 

1  PAUL  GOUGH  ASSOC  INC 

ATTN  P  S  GOUGH 
1048  SOUTH  ST 
PORTSMOUTH  NH  03801-5423 

1  PHYSICS  INTRNTNL  LIBRARY 

ATTN  H  WAYNE  WAMPLER 
PO  BOX  5010 

SAN  LEANDRO  CA  94577-0599 


HERCULES  INC 
ATTN  R  CARTWRIGHT 
AEROSPACE 
100  HOWARD  BLVD 
KENVILLE  NJ  07847 


PRINCETON  CMBSTN  RSRCH  LABS  INC 
ATTN  N  A  MESSINA 
PRINCETON  CORPORATE  PLAZA 
11  DEERPARK  DR  BLDG  IV  SUITE  119 
MONMOUTH  JUNCTION  NJ  08852 


HERCULES  INC 
ATTN  B  M  RIGGLEMAN 
HERCULES  PLAZA 
WILMINGTON  DE  19894 

MARTIN  MARIETTA  ARM  SYS 
ATTN  JIM  TALLEY 
ROOM  1309 
LAKESIDE  AVENUE 
BURLINGTON  VT  05401 

MBR  RESEARCH  INC 
ATTN  DR  MOSHE  BEN  REUVEN 
601  EWING  ST  SUITE  C  22 
PRINCETON  NJ  08540 


3  ROCKWELL  INTRNTNL 
ATTN  BA08 
J  FLANAGAN 
J  GRAY 
R  B  EDELMAN 
ROCKETDYNE  DIV 
6633  CANOGA  AVE 
CANOGA  PARK  CA  91303-2703 

2  ROCKWELL  INTRNTNL  SCIENCE  CTR 
ATTN  DR  S  CHAKRAVARTHY 
DR  S  PALANISWAMY 
1049  CAMINO  DOS  RIOS 
PO  BOX  1085 

THOUSAND  OAKS  CA  91360 


OLIN  CORPORATION 
ATTN  F  E  WOLF 
BADGER  ARMY  AMMO  PLANT 
BARABOOWI  53913 


1  SAIC 

ATTN  M  PALMER 
2109  AIR  PARK  RD 
ALBUQUERQUE  NM  87106 


OLIN  ORDNANCE 
ATTN  E  J  KIRSCHKE 
A  F  GONZALEZ 
D  W  WORTHINGTON 
PO  BOX  222 

ST  MARKS  FL  32355-0222 

ATTN  DR  JOHN  DEUR 
2001  AEROSPACE  PARKWAY 
BROOK  PARK  OH  44142 


1  SOUTHWEST  RSRCH  INSTITUTE 

ATTN  J  P  RIEGEL 
6220  CULEBRA  ROAD 
PO  DRAWER  28510 
SAN  ANTONIO  TX  78228-0510 

1  SVERDRUP  TECHLGY  INC 
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3  TfflOKOL  CORPORATION 
ATTN  R  WELLER 
R  BIDDLE 
TECH  LIBRARY 
ELKTON  DIVISION 
PO  BOX  241 

ELKTON  MD  21921-0241 

1  VERTTAY  TECHLGY  INC 

ATTN  E  FISHER 
A  CRICKENBERGER 
J  BARNES 

4845  MILLERSPORT  HWY 
EAST  AMHERST  NY  14501-0305 

1  UNIVERSAL  PROPULSION  COMPANY 

ATTN  H  J  MCSPADDEN 
25401  NORTH  CENTRAL  AVE 
PHOENIX  AZ  85027-7837 

1  SRI  INTERNATIONAL 

ATTN  TECH  LIBRARY 
PROPULSION  SCIENCES  DIV 
333  RAVENWOOD  AVE 
MENLO  PARK  CA  94025-3493 


ABERDEEN  PROVING  GROUND 


1  CDR  USAATC 

ATTN  STECS  LI  R  HENDRICKSEN 

34  DIR  USARL 

ATTN  AMSRL  WTP  A  HORST 
AMSRL  WT  PA 
T  MINOR 
T  COFFEE 
G  WREN 
ABIRK 
J  DE  SPIRITO 
AJUHASZ 
J  KNAPTON 
CLEVERITT 
M  MCQUAID 
WOBERLE 
P  TRAN 
K  WHITE 
L-M  CHANG 
J  COLBURN 
P  CONROY 
G KELLER 
DKOOKER 
MNUSCA 
T  ROSENBERGER 
AMSRL  WT  PB 
P  PLOSTINS 
M BUNDY 
B  GUIDOS 
D  LYON 
J GARNER 
VOSKAY 
JSAHU 
P  WEINACHT 
HEDGE 
E  FERRY 
AMSRL  WT  PC 
BFORCH 

AMSRL  WT  PD  B  BURNS 
AMSRL  SC  C 
WSTUREK 
C  NIETUBICZ 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes.  Your  comments/answers 
to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  ARL  Report  Number/Author  ARL-TR-1156  (Weinachf) _ Date  of  Report  July  1996 _ 

2.  Date  Report  Received _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest  for  which  the  report 

will  be  used.) _ _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source  of  ideas,  etc.) 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved,  operating  costs 
avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate  changes  to 
organization,  technical  content,  format,  etc.) _ 


Organization 

CURRENT  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 

City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  Current  or  Correct  address  above  and  the 
Old  or  Incorrect  address  below. 


Organization 


OLD  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  tape  closed,  and  mail.) 
(DO  NOT  STAPLE) 


DEPARTMENT  OF  THE  ARMY 


OFFICIAL  BUSINESS 


BUSINESS  REPLY  MAIL 

FIRST  CLASS  PERMIT  NO  0001,APG,MD 


POSTAGE  WILL  BE  PAID  BY  ADDRESSEE 


DIRECTOR 

U.S.  ARMY  RESEARCH  LABORATORY 
ATTN:  AMSRL-WT-PB 

ABERDEEN  PROVING  GROUND,  MD  21005-5066 


NO  POSTAGE 
NECESSARY 
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